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Abstract 

The Hubbard model is a "highly oversimplified model" for electrons in a solid 
which interact with each other through extremely short ranged repulsive (Coulomb) 
interaction. The Hamiltonian of the Hubbard model consists of two pieces; -ffhop 
which describes quantum mechanical hopping of electrons, and H[^t which describes 
nonlinear repulsive interaction. Either -ffhop or i^int alone is easy to analyze, and does 
not favor any specific order. But their sum H = i/hop + -ffint is believed to exhibit 
various nontrivial phenomena including metal-insulator transition, antiferromagnetism, 
ferrimagnetism, ferromagnetism, Tomonaga-Luttinger liquid, and superconductivity. 
It is believed that we can find various interesting "universality classes" of strongly 
interacting electron systems by studying the idealized Hubbard model. 

In the present article we review some mathematically rigorous results on the Hub- 
bard model which shed light on "physics" of this fascinating model. We mainly concen- 
trate on magnetic properties of the model at its ground states. We discuss Lieb-Mattis 
theorem on the absence of ferromagnetism in one dimension, Koma-Tasaki bounds on 
decay of correlations at finite temperatures in two-dimensions, Yamanaka-Oshikawa- 
Affleck theorem on low-lying excitations in one-dimension, Lieb's important theorem 
for half-filled model on a bipartite lattice, Kubo-Kishi bounds on the charge and su- 
perconducting susceptibilities of half-filled models at finite temperatures, and three 
rigorous examples of saturated ferromagnetism due to Nagaoka, Mielke, and Tasaki. 
We have tried to make the article accessible to nonexperts by describing basic defini- 
tions and elementary materials in detail. 
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1 Introduction 

According to the textbook of Ashcroft and Mermin, the Hubbard model is "a highly over- 
simplified model" for strongly interacting electrons in a solid. The Hubbard model is a kind 
of minimum model which takes into account quantum mechanical motion of electrons in a 
solid, and nonlinear repulsive interaction between electrons. There is little doubt that the 
model is too simple to describe actual solids faithfully. 

Nevertheless, the Hubbard model is one of the most important models in theoretical 
physics. In spite of its simple definition, the Hubbard model is believed to exhibit various 
interesting phenomena including metal-insulator transition, antifcrromagnetism, ferrimag- 
netism, ferromagnetism, Tomonaga-Luttingcr liquid, and superconductivity. Serious theo- 
retical studies have also revealed that to understand various properties of the Hubbard model 
is a very difficult problem. We believe that in course of getting deeper understanding of the 
Hubbard model, wc will learn many new physical and mathematical techniques, concepts, 
and ways of thinking. Perhaps a more important point comes from the idea of "universal- 
ity." We believe that nontrivial phenomena and mechanisms found in the idealized Hubbard 
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model can also be found in other systems in the same "universality class" as the ideahzed 
model. The universality class is expected to be large and rich enough so that it contains 
various realistic strongly interacting electron systems with complicated details which are 
ignored in the idealized model. 

The situation is very similar to that of the Ising model for classical spin systems. The 
Ising model is too simple to be a realistic model of magnetic materials, but has turned out to 
be extremely important and useful in developing various notions and techniques in statistical 
physics of many degrees of freedom. Many important universality classes (of spin systems 
and field theories) were discovered by studying the Ising model. 

In the present article, we review some mathematically rigorous results^ known for the 
Hubbard model. We shall concentrate ourselves mainly on magnetic properties of the model 
at its ground state, i.e., at the zero temperature. We have also decided not to cover many 
important rigorous and/or exact results in one-dimensional models based on the Bethe ansatz 
solutions. Even with these restrictions, we do not try to cover all the existing rigorous 
results. We recall that there is an excellent review article by Lieb [1] which covers wider 
topics than we do here. As for the more restricted topics of Nagaoka's ferromagnetism, flat- 
band ferromagnetism, and some related topics, there is a separate review [2] which is more 
detailed and elementary than the present one. 

2 Hubbard Model 

2.1 Definition of the Hubbard Model 

We first give general definition of the Hubbard model^. Let the lattice A be a collection 
of sites x,y,.... Physically speaking, each lattice site corresponds to an atomic site in 
a crystal. In the standard Hubbard model, one simplifies the situation considerably, and 
assumes that each atom has only one electron orbit and the corresponding orbital state is 
non-degenerate^. Of course actual atoms can have more then one orbits (or bands) and 
electrons in the corresponding states. The philosophy behind the model building is that 
those electrons in other states do not play significant roles in low-energy physics that we are 
interested in, and can be "forgotten" for the moment. See Figure 1. 

By cj.^, we denote the operator which creates an electron with spin a =t, i at site x & A. 

^In order to reduce the number of references, we have decided not to include many important references 
on the related topics which do not provide rigorous results. 

^The readers who are new to the filed are recommended to take a look at [2], which contains more careful 
introduction to the Hubbard model. 

^Such a model is usually referred to as a single-band Hubbard model. This terminology is confusing since 
such a model can possess more than one (single-electron) band depending on the lattice structure. Perhaps 
"single-orbital Hubbard model" is a better terminology. 
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Figure 1: A highly schematic figure which explain philosophy of tight- binding descriptions, 
a) A single atom with multiple electrons in different orbits, b) When atoms get together 
to form a solid, electrons in the black orbits become itinerant, while those in the light gray 
orbits are still localized at the original atomic sites. Electrons in the gray orbits are mostly 
localized around the atomic sites, but tunnel to nearby gray orbits with a non-negligible 
probabilities, c) We only consider the electrons in the gray orbits, which are expected to 
play essential roles in determining various low energy physics of the system, d) If the gray 
orbit is non-dcgcnerate, we get a lattice model in which electrons live on lattice sites and 
hop from site to another. 
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The corresponding annihilation operator is Cx^a-, and nx,c — c^x,a^x,(T is the number operator. 
These fermion operators obey the canonical anticommutation relations 

{''sCjO-' '^J/,r} — ^x,ySa,T-) (2-1) 

and 

{<^i,(T) ^]/,t} ~ {Cx,a, Cy^r} = 0, (2-2) 

where {A, B} = AB + BA. 

By $vac we denote the state without any electrons. We have c^^o-'^' vac = for any x e 
A and a =t) i- The Hilbert space of the model is generated by the states obtained by 
successively operating the creation operator cj.^ with various x and a onto the state $vac- 
Since the anticommutation relation (2.2) implies {c\.,^Y — 0, each lattice site can either be 
vacant, occupied by an | or | electron, or occupied by both | and J, electrons. The total 
dimension of the Hilbert space is thus^ 4l^l. 

The Hamiltonian of the Hubbard model is most naturally represented as the sum of two 
terms as 

H — Hhop + Hint- (2.3) 
The most general form of the hopping Hamiltonian i?hop is^ 

Hhop — X! X! ^3:,yCl:,a'^y,a- (2-4) 

x,yeAcT='\,i 

The hopping amplitude t^^y = ty^^, which is assumed to be real, represents the quantum 
mechanical amplitude that an electron hops from site x to y (or from y to x). When x = y, 
the summand in (2.4) becomes tx^xcl:a^x,a = ix,x^x,a) which is nothing but a single-body 
potential. 

The interaction Hamiltonian i^int is written as 

where t/j, > is a constant. The Hamiltonian represents a nonlinear interaction which 
raises the energy by Ux when two electrons occupy a single orbital state at x. Although 

the original Couloinl) iiiteiaetioii is long ranged, we have "oversimplified" the situation and 

^Throughout the present article we denote by l^l the number of elements in a set S. 

^The standard convention is to put a minus sign in front of the summation in (2.4), and to assume 
ix,y > 0. However it seems that there is no simple reason that the hopping amplitude should have such 
signs. If the system is bipartite (See Definition 5.1), one can change the sign of all tx,y (x ^ y) by performing 
a gauge transformation c|. ^ — cj. ^ for all x G A. 
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took into account the strongest part out of the interaction^. Another interpretation is that 
the Coulomb interaction is screened by the electrons in different orbital states which we had 
decided to forget. 

2.2 Some Physical Quantities 

We shall define some basic conserved quantities. The total number operator 

A'e = EK,T + ^a.,i) (2-6) 
xeA 

commutes with the Hamiltonian H. Although there are some conserved quantities other than 
Ne, one usually discusses stationary states or equilibrium states of the system by keeping 
the eigenvalue or the expectation value of A'^e constant*. In the present article, we mostly 
consider^ the Hilbert space in which the number operator has a fixed eigenvalue A'c- Since 
each lattice site can have at most two electrons, we have < < 2|A|. The total electron 
number is the most fundamental parameter in the Hubbard model. 
The spin operator = {S^^\S^^\ S^^) at site x is defined as 

si^'^-l E 4,.(P^"^,.c.,., (2.7) 

for a = 1,2, and 3, where p^^^ are the Pauli matrices. The operators for the total spin of 
the system are defined as 

sti^T.Si'-K (2.8) 

for q; = 1,2, and 3. The operator S^^t commute with both the hopping Hamiltonian i^hop 
(2.4) and with the interaction Hamiltonian Hi^t (2.5). In other words, these Hamiltonians 
are invariant under any global rotation in the spin space. 

As the operators S^^t with a = 1,2,3 do not commute with each other, we follow the 
convention in the theory of angular momenta, and simultaneously diagonalize the total spin 
operators Sl^ol, (Stot)^ = J2a=i{^tot)'^ ^ ^ind the Hamiltonian H. We denote by 5*^^ and 

^There are many important works on various extended Hubbard models in which one takes into account 
other short range interactions which arise from the original Coulomb interaction. See [3, 4, 5, 6] and many 
references therein. 

^For example the total spin is also a conserved quantity. But we do not fix its eigenvalue or expectation 
value, since the total spin is not definitely conserved in reality because there is an LS coupling and the actual 
solids are not rotation invariant. The situation for Ne is essentially different since the charge conservation 
is an exact law. See Section 2.2 of [2]. 

^Section 5.4 is the only exception. 
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'S'tot('S'tot + 1) the eigenvalues of Si^t and (Stot)^, respectively. For a given electron number 
A^e, we let 

_ r NJ2 when < iVe < |A|; . . 

'^'"^ ~ I |A| - (iVe/2) when |A| < iVe < 2|A|. ^ ^ 

Then the possible values of ^tot are ^tot = 0, 1, . . . , S'^ax (or Stot = 1/2, 3/2, . . . , Sr^a^). 

When we discuss magnetism of the system, the most important issue is to determine the 
value of Slot in the ground state(s). If the total spin of the ground state grows proportionally 
to the number of sites |A| as we increase the size of A, we say that the system exhibits 
ferromagnetism in a broad sense. This roughly means that the system behaves as a "magnet." 
If the total spin of the ground state (s) coincides with the maximum possible value ^max, we 
say that the system exhibits saturated ferromagnetism. 

The following quantity will be useful in the later analysis. 

Definition 2.1 (The lowest energy for each S'tot) Fix the electron number N^. For S = 
0, 1, . . . , 5'inax (or S = 1/2, 3/2, . . . , S^ax), we denote by £'min(>S') the lowest possible energy 
among the states which satisfy — N^^ and (Stot)^^ — S{S + 1)$ (i.e., Stot — S). 

The appearance of saturated ferromagnetism is equivalent to have £^min('S') > -Eniin('S'niax) 
for any S such that S < S^aax- 

3 Basic Facts about the Model 

In order to understand the meaning of the Hamiltonian of the Hubbard model, we discuss 
physics we encounter in two limiting situations. 

3.1 Non-Interacting System 

Let us assume that the Coulomb interaction in i^int (2.5) satisfies Ux — for any x e A. Since 
the remaining Hamiltonian H = i^hop (2-4) is a quadratic form in fermion operators, it can 
be diagonalized easily (in principle). The single-electron Schrodinger equation corresponding 
to the hopping Hamiltonian i^hop (2.4) is 

2/eA 

where (p — {(fx)xeA is a single-electron wave function, and e is the single-electron energy 
eigenvalue. We shall denote the eigenvalues and the eigenstates of (3.1) as €j and </5*--'^ = 
iVx '')xeA, respectively, where the index takes values j = 1, 2, . . . , |A|. We count the energy 
levels with taking degeneracies into account, and order them as ej < e^+i. 
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Let us discuss a simple and standard example. Take a one dimensional lattice A = 
{1,2,..., N}, and impose a periodic boundary condition which identifies the site 1 with the 
site + 1. As for the hopping matrix elements, we set t^^x+i = ^rr+i,x = —t, and t^^y = 
otherwise. The corresponding Schrodinger equation (3.1) can be solved easily. By using the 
wave number k = 2'Kn/N (with n = 0, ±1, ±2, . . . , ± {(A^/2) - 1} , A^/2), the eigenstates and 
the eigenvalues can be written as N~^/'^ exp[ikx\ and e{k) — —2tcosk, respectively. If one 
makes a suitable correspondence between n and j — 1,2, . . . , N, we get the desired energy 
level Ej. 

We return to the general setting, and define fermion operators corresponding to the 
eigenstates (^^^^ = ((^(f))a;gA by 

4. = E V'^K^- (3-2) 
xeA 

By using the orthonormality of the set of eigenstates (¥''•"'•' )j=i.2,...,|A| (we redefine the eigen- 
states if they do not form an orthonormal set), one finds that the inverse transformation of 
(3.2) is Cx^a = J2''j^i'Px^cij,cr- Substituting this into (2.4), and by using (3.1), we find that 
-f^hop can be diagonalized as 

Here be interpreted as the electron number operator for the j-th single- 

electron eigenstate. 

Let A, B be two arbitrary subsets of {1, 2, ... , |A|} which satisfy \A\ -\- \B\ — N^. By 
using (3.3), we find that the state 

^A,B= \Y{a]A fllSurvac (3.4) 

is an eigenstate oi H — i^hop and its energy eigenvalue is 

EA,B^j:^j + E^r (3-5) 
jeA jeB 

By choosing subsets A, B which minimize Ea,b, we get ground state (s) of the non- interacting 
model. 

In particular if the corresponding single-electron energy eigenvalues are nondegenerate, 
i.e., Ej < Sj+i, and is even, the ground state of H = i^hop is unique and written as 

/Nc/2 \ 

*GS= I jgajXj^va.. (3.6) 
8 
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Figure 2: Schematic picture of the ground state of a non-interacting many electron system. 
The lowest Ne/2 single-electron energy levels are "filled" by both up spin and down spin 
electrons. The state naturally exhibits paramagnetism known as Pauli paramagnetism. 



This is nothing but the state obtained by "filling up" the low energy levels by up and down 
spin electrons, as one learns in elementary quantum mechanics (Figure 2). It is easily verified 
that the above state has a definite total spin S'tot = 0. The ground state (3.6) exhibits no long 
range order. A system with no magnetic ordering is usually said to exhibit paramagnetism^*^. 

In the simple example in one dimension, all the energy levels except £(0) = —2t and 
e{Tr) = 2t are two-fold degenerate. In this case the ground state oi H = //hop may not 
be unique for some values of N^.. However the degeneracy of the ground states is at most 
four-fold, and the total spin of the ground states can take values in S'tot = 0, 1/2, or 1. We 
can conclude that the property of the ground state (s) is essentially the same as the models 
without degeneracy. In general we can draw the same conclusion unless the single-electron 
spectrum has a bulk degeneracy. 

In a single-electron eigenstate of the example in one dimension, the electron is in a plane 
wave state with a definite wave number k. The fact that the Hamiltonian H = i^hop is 
diagonalized as (3.3) implies that the electrons behave as "waves" in this non-interacting 
(Hubbard) model. The same comment applies to any translation invariant (Hubbard) model 
with Ur = 0. 



3.2 Non-Hopping System 

Let us next assume that the hopping matrix elements in i?hop (2-4) satisfy t^^y — for any 
x,y & A. Then the remaining Hamiltonian H — H-^^t (2-5) is already in a diagonal form. A 

-^•^More precisely, this is true when one talks only about magnetism carried by electron spins. If one takes 
into account magnetism induced by orbital motion of electrons, the system may exhibit diamagnetism. 
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general eigenstate can be written as 

*^.^-fn4,T)fn4,i)*va.. (3.7) 

\xex J \x& / 

Here X and Y are arbitrary subsets of A, and represent lattice sites which are occupied by 
up-spin electrons and down-spin electrons, respectively. The total electron number in this 
state is A^e = l-'^l + 1^1, and the energy eigenvalue is given by 

Exy^ E (3-8) 

The ground state for a given electron number Ag can be constructed by choosing subsets 
X, Y that minimize the energy Ex^y- When one has < |A|, one can always choose X and 
Y such that X (lY — ^. Thus the ground state has energy equal to 0. 

The ground states of the non-hopping Hubbard model has no magnetic long range order. 
Again the system is paramagnetic. It is also clear (from the beginning) that the electrons 
behave as "particles" in non-hopping models. 

3.3 Hubbard Model is Difficult, but it is Interesting 

We have investigated the properties of the two pieces i^hop and ifint in the Hubbard Hamilto- 
nian. It turned out that both i/hop and iifint are easy to analyze. We also found that neither 
of them favor any magnetic ordering. 

We also observed, however, that electrons behave as "waves" in i^hop, while they behave 
as "particles" in i^int. How do they behave in a system with the Hamiltonian which is a sum 
of these totally different Hamiltonians? This is indeed a fascinating problem which is deeply 
related to the wave-particle dualism in quantum physics. We might say that many of the 
important models in many-body problems, including the ip'^ quantum field theory and the 
Kondo problem, are minimum models which take into account at the same time the wave- like 
nature and the particle-like nature (through point-like nonlinear interaction) of matter. 

From a technical point of view, the wave-particle dualism implies that the Hamiltonians 
i^int are -ffhop do not commute with each other. Even when each Hamiltonian is diagonal- 
ized, it is still highly nontrivial (or impossible) to find the properties of their sum. Of course 
mathematical difficulty does not automatically guarantee that the model is worth studying. 
A really exciting thing about the Hubbard model is that, though the Hamiltonians ifhop and 
iJint do not favor any nontrivial order, their sum H — Hhop + Hmi is believed to generate 
various nontrivial order including antiferromagnetism, ferromagnetism, and superconductiv- 
ity. When we sum up the two innocent Hamiltonians i^hop and i^mt, competition between 
wave-like character and particle-like character (or between linearity and nonlinearity) takes 
place, and one gets various interesting "physics." To confirm this fascinating scenario is a 
very challenging problem to theoretical and mathematical physicists. 
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4 Results for Low Dimensional Models 



We discuss some theorems which are proved by using special natures of low-dimensional 
systems. 

4.1 Lieb-Mattis Theorem 

Theorems discussed in the present and the next sections state that the Hubbard model 
does not exhibit interesting long range order under some conditions. The main purpose of 
studying an idealized model like the Hubbard model should be to show that some interesting 
physics do take place. Results which say something does not take place may be regarded 
as less exciting. But to have a definite knowledge that something does not happen under 
some conditions is very useful and important even if our final goal is to show something does 
happen. 

The classical Lieb-Mattis theorem [7] states (among other things) that one can never have 
ferromagnetism in one-dimensional Hubbard model with only nearest neighbor hoppings^^. 

Theorem 4.1 (Lieb-Mattis theorem) Consider a Hubbard model on a one- dimensional 

lattice A = {1, 2, ... , N} with open boundary conditions. We assume that the hopping matrix 
elements satisfy \tx,y\ < oo when x = y , < \tx,y\ < oo when \x — y\ = 1, and are vanishing 
otherwise. We also assume \Ux\ < oo for any x E A. Then the quantity £^mm(>S') (see 
Definition 2.1) satisfies the inequality 

E^US)<E^US + l), (4.1) 
for any S = 0,1, ... , S'max - 1 (or S = 1/2,3/2, . . . , S'max - I)- 

One of the most important consequences of the Lieb-Mattis theorem is that any one- 
dimensional Hubbard model in the above class has the total spin ^tot = (or S'tot = 1/2) 
in its ground state. One cannot conclude from this fact alone that the system exhibits 
paramagnetism, but can conclude that there is no ferromagnetism. 

Theorem 4.1 docs not apply to models with periodic boundary conditions. But it seems 
reasonable that boundary conditions do not change the essential physics provided that the 
system is sufficiently large. If there exist hoppings to sites further than the nearest neighbor, 
on the other hand, the story is totally different. We do not only find that the proof of 
Theorem 4.1 fails, but we also find essentially new physics. See Section 6.5. 

Theorem 4.1 is proved by noting that, in a suitable basis, the Hamiltonian is written as 
a matrix whose nondiagonal elements are nonpositive, and by using the standard Perron- 
Frobenius argument. Similar argument was used by Lieb and Mattis in their study of the 
Heisenberg quantum spin system [8]. 

^^The present theorem appears in the Appendix of [7]. The main body of [7] treats interacting electron 
systems in continuous spaces. 
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4.2 Decay of Correlations at Finite Temperatures 

Among other rigorous results which show the absence of order is the extension by Ghosh [9] 
of the wellknown theorem of Mermin and Wagner. Ghosh proved that the Hubbard model in 
one- or two-dimensions does not exhibit symmetry breaking related to magnetic long range 
order at any finite temperatures. By using the same method, one can also prove the absence 
of superconducting U{1) symmetry breaking. 

Koma and Tasaki proved essentially the same facts in terms of explicit upper bounds for 
various correlation functions [10]. Among the results of [10] is the following. 

Theorem 4.2 (Koma- Tasaki bounds for correlations) Consider an arbitrary Hubbard 
model in one- or two- dimensions with finite range hoppings. Then there are constants a, 7, 
and we have 



for sufficiently large \x — y\, where {■ ■ ■) ^ denotes the canonical average in the thermodynamic 
limit at the inverse temperature (5. Here f{(3) is a decreasing function of P and behaves as 
f{(5) ~ 1//3 for /3 » 5 and f{(5) ^ {2/5)\lnf3\ for (5 ^^5, where 5 is a constant. 

The bounds (4.2) and (4.3) establishes the widely accepted fact that there can be 
no superconducting^^ or magnetic long-range order at finite temperatures in one- or two- 
dimensions. The method employed in [10], i.e., a combination of the McBryant-Spencer 
method and the quantum mechanical global U{1) gauge invariance, is rather interesting. It 
is amusing that only by using the U{1) symmetry which exists in any quantum mechanical 
system, one gets upper bounds for correlations which are almost optimal at low temperatures 
(especially in one-dimension). 

4.3 Yamanaka-Oshikawa- Affleck Theorem 

We discuss a recent important theorem by Yamanaka, Oshikawa, and Affleck [11, 12] about 
low-lying excitations in general electron systems on a one-dimensional lattice. The theorem 
is an extension of Lieb-Schultz-Mattis theorem for quantum spin chains. It can also be in- 
terpreted as a nonperturbative version of Luttinger's "theorem" restricted to one-dimension. 

^^One can easily extend (4.2) to rule out the condensation of other types of electron pairs. 




(4.2) 



and 




(4.3) 
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We consider the Hubbard modeP^ on the one-dimensional lattice A = {1, 2, ... , N} with 
periodic boundary conditions. The model is characterized by positive integers R and P, 
which are the range of hopping and the period of the system, respectively. We assume 
tx,y = whenever |a; — 7/| > R, and tx+p,y+p = tx,y, Ux+p = Ux for any x,y E A. Under this 
general assumption, we have the following. 

Theorem 4.3 (Yamanaka-Oshikawa- Affleck theorem) Consider the infinite volume limit 
iV ^ oo with a fixed electron density u = N^/N. If Pu/2 is not an integer, then we have 
one of the following two possibilities. 

i) There is a symmetry breaking, and the infinite volume ground states are not unique. 

ii) There is a gapless excitation above the infinite volume ground state. 

In other words, the theorem rules out the third possibility; 
Hi) The infinite volume ground state is unique, and there is a finite gap above it. 

Note that iii) does take place if Pu/2 is an integer and the system describes an innocent 
insulator In a non-interacting system, it is evident that iii) is impossible when Pu/2 is 
not an integer, since there is a partially filled band. The above theorem guarantees that we 
cannot change the situation by introducing strong interaction. This is of course far from 
obvious. 

The proof of Theorem 4.3 is based on the following exphcit construction of a (trial) low- 
lying excitation^^ . Consider a system on a periodic chain of length A^, and assume that the 
ground state <I>gs is unique. We define the "twist" operator hy U = eiqi[2niYl,x=i{x / N)nx^{\ 
which introduces a gradual twist in the U{1) phase of the up-spin electron field, and consider 
the trial excited state vl* = U^^qs- It is not hard to show that - Eqs = 0{l/L). 

Thus ^ contains a low-lying excited state provided it is orthogonal to the ground state. To 
see the orthogonality, let Tp be the translation by P, and assume that the ground state is 
chosen as Tp$Gs = '^'gs- It is easily found that Tp^ = e'^'"^^'^ ., and hence the trial state ^ 
is orthogonal to $gs whenever Pu/2 is not an integer. 

The above result also implies that, in the case ii), the gapless excitation has a crystal 
momentum k^, — ttu. If we interpret this excitation as obtained from the ground state by 
moving an electron at a "fermi point^^" to the other fermi point, we find that the fermi 
momentum kp satisfies = 2k-F, and hence kp = tiv/2. This is nothing but the fermi 
momentum of the free system. Therefore the fermi momentum is not "renormalized" by 
strong interaction among electrons. As far as we know, this is the only general rigorous 

^•^The theorem applies to a much larger class of lattice electron systems. It is especially meaningful when 
applied to the Kondo lattice model. 

^''Consider, for example, a two-band system with a band gap which has P = 2. Then, in a free system, 
1/ = 1 (the half-filling) corresponds to an insulator with a charge gap. 

^^An extra care is needed to discuss infinite volume limits [12]. 

^^In a one- dimensional model, fermi surface (if any) becomes two "fermi points." 



13 



result which gives precise meaning to fermi momentum in truly interacting many electron 
systems. 



5 Half-Filled Systems 



A system in which the electron number is identical to the number of sites |A| is said to 
be half-filled, since the maximum possible value of A^e is 2|A|. The system becomes half- 
filled if each atom provides one electron to the system. Thus the half-filled models represent 
physically natural situations. Half-filled models have nice properties'^ from mathematical 
point of view as well, and there are some very nice rigorous results. 

5.1 Perturbation for U ^ t 

Let us first look at the ground states of the non-hopping model with i^hop — 0. Wc here 
assume that Ux > for any a; G A. As we found in Section 3.2, one can choose X fl F = 
in the state to get a ground state with Ex,y = 0, provided that < |A|. Since we 

have Ae = |A|, the assumption X (lY — ^ automatically imphes X L)Y — A. Therefore the 
ground state ^x,y (3.7) with X r\Y — $ can be written also as 



where cr — {cr{x))x&A is a collection of spin indices a{x) J,. By using the terminology of 
spin systems, one can call cr a spin configuration. One can say that the degeneracy of the 
ground states (5.1) precisely corresponds to all the possible spin configurations. 

Let us take into account effects of nonvanishing ifhop by using a simple-minded perturba- 
tion theory. As the diagonal elements of ifhop, i-e., T,x,a tx,x cl^^c^^a = T,x tx,x{^x,^+nx,i), only 
shift the energy of the states by a constant amount (independent of «t) , they can be omit- 
ted in the lowest order perturbation calculation. Let us denote by Hhop = J2xjty,(7 ^x,y cl:^a^y,a 
the off-diagonal part of if hop- By operating ifhop once onto \I'cr, an electron moves, and we 
get a state with one vacant site and one doubly occupied site. The resulting state is not a 
ground state of Hi^i. We thus find that the lowest order contribution from this perturbation 
theory comes from the second order. 

Figure 3 shows a process that is taken into account in the second order perturbation 
theory. The electron at site x hops to site y with the transition amplitude tx,y, and generates 
a new state with extra energy Uy. Then one of the two electrons at site y will hop back to 

^''Some half-filled models can be mapped onto a Hubbard model with attractive interaction via a partial 
hole-particle transformation. This fact plays crucial role in the proof of Lieb's theorem. 




(5.1) 
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X y X y X y 

Figure 3: When electrons hop twice, spins on sites x and y may be exchanged. This is the 
ultimate origin of antiferromagnetic nature in the half-filled Hubbard model. 

site X, and we recover one of the ground states. In this process, spins at the sites x and y 
may be exchanged as Figure 3 shows. The hopping between sites x and y is inhibited by the 
Pauli principle if the electronic spins on these two sites are pointing in the same direction. 
We find that this second order perturbation process lowers the energy of states in which the 
spins at sites x and y are not pointing in the same direction (or more precisely, the states in 
which the total spin is vanishing). 

Let us rederive this result in a more formal manner. Let Pq be the projection operator 
onto the subspace spanned by the states \l/cr (5.1) with all the possible <t. That the first 
order perturbation has no contribution can be read off from the fact that Po-f^hop-Po = 0. 
To find out how the degeneracy in the (unperturbed) ground states (5.1) is lifted can be 
determined by the effective Hamiltonian 

HeS = — -Po-f^hopTT — -f^hop-Po 
-Hint 

= -Po S X! "^a;,!/ (^x '^y ~ j)\ -^0- (5.2) 

Here the exchange interaction parameter is given by J^^y = {itx,y)'^ /Ux} + {(tx,yy /Uy}. 
Note that (5.2) is nothing but the Hamiltonian of the S" = 1/2 antiferromagnetic quantum 
Heisenberg spin system. This suggests that the low-energy behavior of the half-filled Hubbard 
model is well described by the Heisenberg antiferromagnets when Ux are much larger than 

tx,y 

5.2 Lieb's Theorem 

In 1989, Lieb proved an important and fundamental theorem for the half-filled Hubbard 
model. The theorem provides, among other things, partial support to the conjecture about 
the similarity of the half-filled Hubbard model and the Heisenberg antiferromagnets. Let us 
first introduce the notion of bipartiteness. 

Definition 5.1 (Bipartiteness) Consider a Hubbard model ( or other tight-binding electron 
model) on a lattice A with hopping matrix elements {tx,y)x ^a- The system is said to be 
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bipartite if the lattice A can be decomposed into a disjoint union of two sublattices as A — 
AU B (with An B = (/)), and it holds that t^^y — whenever x,y&Aorx,yEB. In other 
words, only hoppings between different sublattices are allowed. 

Then Lieb's theorem [13] for the repulsive Hubbard model is as follows. 

Theorem 5.2 (Lieb's theorem) Consider a bipartite Hubbard model. We assume |A| is 
even, and the whole A is connected^^ through nonvanishing t^^y. We also assume Ux — U > 
for any x & A. Then the ground states of the model are nondegenerate apart from the trivial 
spin degeneracy^^ , and have total spin Stot = ll^l ~ 1-^11/2- 

The total spin S'tot of the ground state determined in the theorem is exactly the same 
as that of the ground state (s) of the corresponding Heisenberg antiferromagnet on the same 
lattice. In fact the conclusion of the theorem is quite similar to that of the Lieb-Mattis 
theorem [8] for Heisenberg antiferromagnets. However the straightforward Perron-Probenius 
argument used in the proof of the latter theorem does not apply to the Hubbard model 
except in one dimension. (See 4.) This is not only a technical difficulty, but is a consequence 
of the important fact that quantum mechanical processes allowed in the Hubbard model 
are in general much richer and more complex than those in the Heisenberg model. Lieb's 
proof is compactly presented in a Letter, but is deep and elegant. The proof again makes 
use of a kind of Perron-Frobcnius argument, but is based on an interesting technique called 
spin-space reflection positivity. 

Lieb's theorem is valid for any value of Coulomb repulsion U, only provided that it is 
positive. It is quite likely that physical properties of the Hubbard model are drastically 
different in the weak coupling region with small U and in the strong coupling region with 
large U. It is very surprising and interesting that a single proof of Lieb's covers the whole 
range with U > and clarifies the basic properties of the ground states. 

It should be noted, however, that the knowledge of the total spin of the ground states 
in finite volume does not necessarily determine the properties of the ground states in the 
corresponding infinite system. When two sublattices have the same number of sites as 
\A\ = I S I, for example, one knows that finite volume ground state is unique and has (S'tot — 0. 
Although one might well conclude that the system has no long range order in its ground 
states, this is not true. It is certainly possible that infinite volume ground states exhibit 
long range order and symmetry breaking (such as antiferromagnetism of superconductivity) 
even when finite volume ground state is unique and symmetric. (See, for example, [14].) 

^^More precisely, for any x,y G A, one can find a sequence of sites a;o, xi, . . ., xn with xq = x, xn = y, 
and 7^ for i = 0, 1, . . . , TV - 1. 

^^In any quantum mechanical system with a rotation invariant Hamiltonian, an eigenstate of the Hamil- 
tonian with the angular momentum J is always (2 J + l)-fold degenerate. 
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By knowing any finite volume ground state has S'tot = 0, however, one can rule out the 
possibility of ferromagnetism. 

By using Lieb's results in [13], one gets some information about excited states. For 
example, by combining Theorem 1 in [13] with the method of [8], one can easily prove the 
inequahty^° 

^min(^)<£^min(^+l) (5.3) 

for any \\A\ - \B\\/2 <S< (|A|/2) - 1. 

Another theorem which suggests the similarity between the half-filled Hubbard model and 
the Heisenberg antiferromagnets is the following, proved by Shen, Qiu, and Tian [15, 16] by 
extending Lieb's method. 

Theorem 5.3 (Explicit signs of correlations) Assume the conditions for Theorem 5.2. 
If we denote by $gs the ground state of the model, we have the inequalities 

<b Q Q fl) when x,y e A, or x,y e B; , . 

*GS,S.-S,$cs;|^Q ^henxeA,yeB,orxeB,yeA, ^^'^^ 

were (■, ■) denotes the quantum mechanical inner product. 

We see that spins on different sublattices have negative correlations, indicating antifer- 
romagnetic tendency. 

The S — 1/2 Heisenberg antiferromagnet on the cubic lattice, for example, is proved 
to exhibit an antiferromagnetic long range order at sufficiently low temperatures or in the 
ground states [17, 18]. It is likely that the same statements hold for the half- filled Hubbard 
model with sufficiently large U. But, for the moment, there are no methods or ideas which 
are useful in proving the conjecture. To extend the powerful (but not very natural) method 
of [17, 18] based on the (spatial) reflection positivity seems hopeless. 

In the Hubbard model and related models at half-fiUing, there have been proved several 
interesting general results. Among the recent examples arc the uniform density theorem [19] , 
the solution of the flux phase problem [20, 21], and the stability of the Peierls instability 



5.3 Lieb's Ferrimagnetism 

A very important corollary of Lieb's theorem 5.2 is that the half-filled Hubbard models 
on asymmetric bipartite lattices universally exhibit a kind of ferromagnetism (in the broad 
sense), or more precisely, ferrimagnetism [13]. 

Take, for example, the so called CuO lattice in Figure 4. The lattice can be decomposed 
into two sublattices distinguished by black sites and white sites. When the black sites form 

wish to thank Shun-Qing Shen and ElUott Lieb for discussions related to this corollary. 
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Figure 4: An example (the so called CuO lattice) of a bipartite lattice in which the number 
of sites in two sublattices are different. Lieb's theorem implies that the half-filled Hubbard 
model defined on this lattice exhibits ferrimagnetism. 

a square lattice with side L, there are black sites and white sites. We define the 
Hubbard model on this lattice, and put nonvanishing hopping t^^y on each bond in the 
lattice, and put Coulomb interaction t/ > on each site. Then Lieb's theorem implies that 
the ground state of this Hubbard model has total spin S'tot = ||^| — |-B||/2 = Since 
the total spin magnetic moment of the system is proportional to the number of lattice sites 
3L^, we conclude that the model exhibits ferromagnetism in the broad sense. 

Of course the present ferromagnetism is not a saturated ferromagnetism in which all the 
spins in the system completely align with each other. As the inequality (5.4) suggests, spins 
on neighboring sites have tendency to point in the opposite direction. But the big difference 
in the numbers of sites in the sublattices cause the system to possess bulk magnetic moment. 
Such a magnetic ordering is usually called ferrimagnetism^^. 

One can similarly construct Hubbard models which exhibit ferrimagnetism on any bipar- 
tite lattice in which the difference in the number of sites in two sublattices is proportional 
to the system size. The value of ?7 > is again arbitrary, so Lieb's ferrimagnetism covers 
surprisingly general class of models including weakly coupled ones as well as strongly coupled 
ones. 

If one recalls the conclusion of Section 3.1 that systems with U = Q exhibit paramag- 
netism, one might feel it somehow contradicting that the above ferrimagnetism appears for 
arbitrarily small [/ > 0. This is one of the special features of Lieb's ferrimagnetism. In 
the single-electron Schrodinger equation corresponding the Hubbard model on Figure 4, for 
example, the eigenstates for the eigenvalue e = are L^-fold degenerate. (The eigenvalue 

^^It is also possible to consider order parameters to see the order is indeed ferrimagnetic [15]. 
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£ = is at the middle of the single-electron spectrum.) Consequently the ground states of 
the half-filled (A^e — \M — 3L^) system with U = are highly degenerate, and the total spin 
can take values in Stot = 0, 1, . . . ,L'^/2. The role of Coulomb interaction U is to lift this 
degeneracy, and select states with the largest magnetic moment as ground states. 



5.4 Kubo-Kishi Bounds on Susceptibilities at Finite Tempera- 
tures 

A theorem which can be regarded as a finite temperature version of Lieb's theorem was 
proved by Kubo and Kishi [23]. It deals with the charge susceptibility and the on-site 
pairing (superconducting) susceptibility in a half-filled system at finite temperatures. 

We define the thermodynamic function^^ J corresponding to the grand canonical ensemble 

by 



J(P,l^, {'yx)xeA, {Vx)x&a) 

( 



p 



log Tr exp 



-P 



\ 



, (5.5) 



where /5 and \i are the inverse temperature and the chemical potential, respectively, and the 
trace is taken over the Hilbert spaces with all the possible electron numbers. We added to 
the Hamiltonian two fictitious external fields (7a;)a;eA and {r]x)x^K to test for the possible 
charge ordering and superconducting ordering, respectively. 

We define the charge susceptibility x'^ and the on-site pairing susceptibility by 



and 



d d 

Xq(/^>/^) = - —^^JiP^I^^ (7a;), {Vx)) 



>o. 



(7x)=(J7x)=0 



The Fourier transformation of the external fields are 



> 0. 



(5.6) 



(5.7) 



(7x)=(t?x)=0 



7q 



xeA 



iAr^/'E7.e''^-, ^q = iAr^/^E^.^e^"-^ 

xeA 



(5.8) 



where q is a wave vector corresponding to the lattice A (which we assume to have a periodic 
structure) . 

Then the Kubo-Kishi theorem states that 



22We have J = -pV = F - G. 
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Theorem 5.4 Consider any bipartite (see Definition 5.1) Hubbard model with Ux — U > 
for any x e A. Then for any (3 > and for any wave vector q, we have 

xl{(3,U/2)<^, and x^l/S, ^/2) < |. (5.9) 

Note that the choice ji — U/2 corresponds to the half-fiUing. The theorem states that 
the charge and the on-site paring susceptibihties for any wave vector q are finite in a half- 
filled model at finite temperatures. This means that the model does not exhibit any CDW 
ordering or superconducting ordering. 



6 Ferromagnetism 

Ferromagnetism, where almost all of the spins in the system align in the same direction, is 
a remarkable phenomenon. Standard theories about the origin of ferromagnetism have been 
the Heisenberg's exchange interaction picture, and the Stoner criterion derived from the 
Hartrcc-Fock approximation for band electrons. But there have been serious doubts if these 
theories really explain the appearance of ferromagnetism in a system of electrons interacting 
via spin-independent Coulomb interaction. One of the motivations to study the Hubbard 
model was to understand the origin of ferromagnetism in an idealized situation. 

As we have seen in the previous section, half-filled models have tendency towards antifer- 
romagnetism. In this section we shall concentrate on systems in which the electron numbers 
deviate from half-filling. 



6.1 Instability of Ferromagnetism 

To see that ferromagnetism is indeed a delicate phenomenon, we discuss two elementary 
results which show that the Hubbard model with certain conditions does not exhibit ferro- 
magnetism^^ . 

The following theorem states that there can be no saturated ferromagnetism if the 
Coulomb interaction U is too small in a system with a "healthy" single-electron spectrum. 

Theorem 6.1 (Impossibility of ferromagnetism for small U) Let {£j}j=i,...,Ar denote 
the single- electron energy eigenvalues with Sj < Cj+i as in Section 3.1. IfO<U< e^^ — e\, 
we have?'^ 

-E'min('5'max ~ l) < -E'min ('5'max) ■ (6-1) 

Thus the ground state of the model does not have Stot — S^ax- 

^^Detailed proofs of the results in the present section can be found in [2]. 
^^Note that the fermi energy sn^, — ei is an intensive quantity. 
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When the density of electrons is very low, it is expected that the chance of electrons to 
collide with each other becomes very small. It is likely that the model is close to an ideal 
gas no matter how strong the interaction is, and there is no ferromagnetism. 

This naive guess is easily justified for "healthy" models in dimensions three (or higher). 
The dimensionality of the lattice is taken into account by assuming that there are positive 
constants c, uq, and d, and the single electron energy levels satisfy 



for any n such that n/\A\ < uq. Note that the right-hand side represents the n dependence 
of energy levels in an usual rf-dimensional quantum mechanical system. Then we have the 
following theorem due to Fieri, Daul, Baeriswyl, Dzierzawa, and Fazekas [24]. 

Theorem 6.2 (Impossibility of ferromagnetism at low densities) Suppose we have ifhop 
satisfying (6.2) with positive c, vq, and d > 2. Then there exists a constant Vi > 0, and the 
same conclusion as Theorem 6.1 holds for any U >0 if Ne/\A\ < ui holds. 

That we have a restriction on dimensionality in Theorem 6.2 is not merely technical. In 
a one-dimensional system, moving electrons must eventually collide with each other for an 

obvious geometric reason. Thus a one-dimensional model cannot be regarded as close to 
ideal no matter how low the electron density is. We do not know whether the inapplicability 
of the theorem to o? = 2 systems is physically meaningful or not. 



As a starting point of our study of ferromagnetism, we consider a toy model with two 
electrons on a small lattice. Interestingly enough, some essential features of ferromagnetism 
found in many-electron systems (that we will discuss later in this section) are already present 
in the toy model. 

The smallest possible model in which we can discuss electron interaction and which is 
away from half-filling is that with two electrons on a lattice with three sites. Consider the 
lattice A = {1,2,3}, and put one electron with a =1 and one with a =j,. The hopping 
matrix is defined by ti^2 — ^2,3 = t', and ^1,3 — t. Note that there are two kinds of hoppings 
t and t'. Since the sign of t' can be changed by the gauge transformation C2,a ~C2,a, we 
shall fix t' > 0. Figure 5 shows the lattice and the hopping. For simplicity, we assume there 
is only one kind of interaction, and set Ui = U2 = U3 = U > 0. We have S'max = 1 because 
Ag = 2. Therefore we can say that there appears saturated ferromagnetism if the ground 
state has Stot — 1; i-e., if it is a part of a spin-triplet. 




(6.2) 



6.2 Toy Model with Two Electrons 
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1 2 3 

Figure 5: The lattice and the hopping of the toy model. By considering the system with two 
electrons on this lattice, we can observe some very important aspects of ferromagnetism in 
the Hubbard model. 




Figure 6: Allowed states and transition amplitudes in the toy model with U = oo. The total 
spin of the ground states can be easily read off from this diagram. 
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Let us take the limit U — > oo, in which the effect of interaction becomes most drastic, 
and consider only those states with finite energies. This is equivalent to consider only states 
in which two electrons never occupy a same site. There are six states which satisfy the 
constraint, and they can be written as ^x,y = cl^|C^ j$vac where x,y = 1,2,3, and x ^ y. 
Transition amplitudes between these states are shown in Figure 6. We find that the problem 
is equivalent to that of a quantum mechanical particle hopping around on a ring consisting 
of six sites. The basic structure of the ground state can be determined by the standard 
Perron- Frobenius sign convention^^ . The ground state for t < is written as 

= ^1,2 + $3,2 - a(t, t')^s,i + $2,1 + $2,3 - a(t, t')$l,3, (6-3) 

and that for t > as 

$gs = $1,2 - $3,2 + Pit, i')$3,l - $2,1 + $2,3 - P{t, t')<^l,3, (6.4) 

where a{t,t') and /3{t,t') are positive functions of t and t'. 

To find the total spin of these states, it suffices to concentrate on two lattice sites, say 

sites 1 and 2, and note that = $1,2 + $2,1 H , and $[^3°^ = $1,2 - $2.1 H • It 

immediately follows that^^ $gs°'' has 5'tot = 0, and $gs°'' has 5'tot = 1- A ferromagnetic 
coupling is generated when t > 0! 

Let us look at the mechanism which generates the ferromagnetism. The states $1^2 and 
$2,1 can be found in the upper left and and the lower right, respectively, in the diagram of 
Figure 6. By starting from $1^2 and following the possible transitions, one reaches the state 
$2,1- In other words, electrons hop around in the lattice, and the spins on sites 1 and 2 are 
"exchanged." When t > 0, the quantum mechanical amplitude associated with the exchange 
process generates the superposition of the two states which precisely yields ferromagnetism. 

Let us briefiy look at the cases with finite U. In Figure 7, we plotted -E'min(O) and £'inin(l) 
for the toy model with t = t'/2 as functions of U. (See Definition 2.1.) As is suggested by 
the result in the U ^ 00 limit, we have ferromagnetism in the sense that i?min(0) > -Emin(l) 
when U is sufficiently large. A level crossing takes place at finite U, and the system is 
no longer ferromagnetic for small U. Even in the simplest toy model, ferromagnetism is a 
"nonperturbative" phenomenon which takes place only when U is sufficiently large. 

The only exception is the case with t = t'. See Figure 8. For this parameter value, the 
ground states are degenerate in spin when U — 0. Ferromagnetic state is the only ground 
state for U > 0. 

^^li the transition amplitude between two states is negative {resp., positive), one superposes the two states 
with the same {resp., opposite) signs. 

^^A quick way to find the total spin of the state ^1.2 + $2,1 is to rewrite the state in the "spin language" 

as $1,2 + $2,1 = c| .f4_^$vac + 4,T4,i*vac = cl,T4,i*vac " cl,i4,T*vac = |T)ili)2 " li)ilT>2' ^nd use the 
standard knowledge about the addition of angular momenta. One can easily convince oneself that the state 
has S'tot = 0. 
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Figure 7: The U dependence of £'min(0) (gray curve) and -E'min(l) (black line) in the toy 
model with t = t'/2. We have ferromagnetism in the sense that £^111111(0) > E^[^{1) when U 
is sufficiently large. We find that ferromagnetism is a "nonperturbative" phenomenon. 



E/f 




Figure 8: The U dependence of i?min(0) (gray curve) and -Emin(l) (black line) in the toy model 
with t = t' . Only for this special parameter, we have ferromagnetism E^^-^{Q) > i?min(l) for 
any value oi U > 0. One can regard this case as the simplest example of the fiat-band 
ferromagnetism that we will discuss in Section 6.4. 
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Figure 9: Schematic picture about the origin of Nagaoka's ferromagnetism. When the hole 
hops around the lattice, the spin configuration is changed. For a model with t^^y > 0, the 
hole motion produces a precise linear combination of various spin configurations which leads 
to ferromagnetic state. 

From Figures 7 and 8, we find that the energy i?jnin(l) of ferromagnetic states is indepen- 
dent of U . As we see in the following, this is a general property of ferromagnetic eigenstates 
in the Hubbard model. An arbitrary state ^ which has total spin 5'tot = >S'max can be written 
as a superposition of states which are obtained by rotating the state ^ which consists only 
of t spin electrons. If one operates the interaction Hamiltonian i/int (2.5) onto the state \E', 
one has i^int^ = because rix^i^ ~ ^- Since the interaction Hamiltonian (2.5) is invariant 
under rotation in spin space, we have shown that i/int^ = 0. One might say that states with 
saturated magnetization do not feel Hubbard type interaction at all. This is one of conve- 
nient (but "oversimplified" ) features in discussing saturated ferromagnetism in the Hubbard 
model. 

6.3 Nagaoka's Ferromagnetism 

The transitions between the states in Figure 6 are generated by hoppings of electrons. One 
can also regard that the transitions are caused by hoppings of a single hole^^, which is 
the site without electrons. At least in the limit ?7 — > cxo, one can say that the origin of 
ferromagnetism in the toy model is the motion of a single hole, which mix up various spin 
configurations with proper signs. 

As Nagaoka [25] demonstrated rigorously, there is a class of many electron models in 
which saturated ferromagnetism is generated by exactly the same mechanism. See Figure 9. 

^'^This is different from tlie notion of liole in tlie usual band tlieory. 
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Nagaoka's theorem (in the extended form of [26] whose complete proof can be found in [2]) 
is as follows. 

Theorem 6.3 (Nagaoka's ferromagnetism) Take an arbitrary finite lattice A, and as- 
sume that^^ tx,y > for any x ^ y, and U^ — oo for any x e A. We fix the electron number 

as Nf. = \K.\ — 1. Then among the ground states of the model, there exist states with total 
spin Slot = Sma.x{= Nc/2). If the system further satisfies the connectivity condition, then the 
ground states have Stot = 'S'max(= Ne/2) and are nondegenerate apart from the trivial spin 
degeneracy. 

The connectivity condition is a simple condition which holds on most of the lattices in 
two or higher dimensions, including the square lattice, the triangular lattice, or the cubic 
lattice. To be precise the condition requires that "by starting from any electron configuration 
on the lattice and by moving around the hole along nonvanishing tx,y, one can get any other 
electron configuration." 

Thouless also reached a similar conclusion [27], but Nagaoka's treatment covers a larger 
class of models including non-bipartite systems. The proof of Nagaoka's theorem (especially 
the recent proof in [26]) is surprisingly simple. It essentially uses the Perron- Frobenius 
argument exactly the same as that we used in Section 6.2 to determine the total spin of the 
ground state of the toy model. 

The requirements that U should be infinitely large and there should be exactly one hole 
are indeed rather pathological. Nevertheless, the theorem is very important since it showed 
for the first time in a rigorous manner that quantum mechanical motion of electrons and 
strong Coulomb repulsion can generate ferromagnetism. The conclusion that the system 
which has one less electron than the half-filled model exhibits ferromagnetism is indeed 
surprising. This is a very nice example which demonstrates that strongly interacting electron 
systems produce very rich physics. 

It is desirable to extend Nagaoka's ferromagnetism to systems with a finite U and with 
a finite density of holes. Although more than thirty years have passed since Nagaoka's and 
Thouless's papers, it is still not known if such extensions are possible. There are, however, 
considerable amount of rigorous works which establish that saturated ferromagnetism does 
not take place in certain situations. See, for example, [28, 29, 30, 31, 32, 33]. 

6.4 Mielke's Ferromagnetism and Flat-Band Ferromagnetism 

Let us once again look at the toy model of section 6.2. As is shown in Figure 8, the 
ground state of the model exhibits ferromagnetism for any U > for the special choice of 

^*As we noted in Section 2.1, this sign of tx,y is opposite from the "standard" choice. In bipartite systems 
(such as those on the square lattice or the cubic lattice with nearest neighbor hoppings), one can change the 
sign of tx,y by a gauge transformation. 
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Figure 10: The Hubbard model on the kagome lattice is a typical example which exhibits 
flat-band ferromagnetism. "Kagome" is a Japanese word for a pattern of woven bamboo of 
baskets. 

the parameters t = f > 0. For this choice of parameters, the energy eigenvalues of the 
corresponding Schrodinger equation are ei = 62 = —t', and £3 = 2t'. The single-electron 
ground states are doubly degenerate. As a consequence, the ground states of the two electron 
system with U = are also degenerate, and can have Stot = or ^tot = 1- The degeneracy 
is lifted for f/ > 0, and the ferromagnetic state is "selected" as the true and unique ground 
state. It is crucial here that the dimension of the degeneracy in the single-electron ground 
states (which is two) is the same as the electron number = 2. 

Mielke [34] showed that there is a class of Hubbard models with many electrons which 
show saturated ferromagnetism through a somewhat similar mechanism. Take, for example, 
the kagome lattice of Figure 10, and define a Hubbard model on it by setting t^^y = t > for 
neighboring sites x, y, t^-^y = for other situations, and = U > for any a; G A. It is worth 
mentioning that the kagome lattice of Figure 10 can be regarded as constructed by putting 
together many copies of the lattice used in the toy model (Figure 5). The energy eigenvalues 
of the corresponding Schrodinger equation can be shown to satisfy ei = 62 = ■ ■ ■ = Sm = — 2t, 
and Ej > —2t for j > M. Here the dimension M of the degeneracy of the single-electron 
ground states is given by M = (|A|/3) + 1, and is proportional to the lattice size. 

We shall fix the electron number as = M, the same as the dimension of the degeneracy. 

Let us consider the case with U = first. Let A and B be arbitrary subsets of 
{1,2, . . . ,Ne} which satisfy \A\ + \B\ = Ne, and consider the state (3.4) obtained by 

creating the corresponding single-electron eigenstates. In the present model on the kagome 
lattice, the fermion operator a^ ^. (3.2) creates one of the single-electron ground states with 
the energy e = —2t. This means that we have ifhop'^'A,^ = ~2tAe$yi.B for arbitrary choice 
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of A and B, and hence ^a,b is a ground state oi H — i?hop- We find that the ground states 
are highly degenerate, and can have S'tot = 0, 1, . . . , S'max = M/2 (or S'tot = 1/2 - • • • j •S'max)- 

What is the effect of nonvanishing Coulomb interaction U in such a situation? Let us 
denote by $| the state obtained by setting A = {1, 2, ... , A^e} and 5 = in $a,_b- Of course 
$1 is one of the ground states of i/hop- As we discussed at the end of Section 6.2, the state $| 
which consists only of t spin electrons "does not feel" Hubbard type Coulomb interaction. 
This means that we have i^int^f = 0. Since is the minimum possible eigenvalue of i^int, 
we find that the state is a ground state of the total Hamiltoniani? — H^op + Hmt for any 
t/ > 0. 

These are all simple observations. A really interesting problem is whether there can be 
ground states other than when [/ > 0. The following theorem due to Mielke shows that 
the ferromagnetic state is indeed "selected" as the true ground state exactly as in the toy 
model. 

Theorem 6.4 (Mielke's flat-band ferromagnetism) Consider the Hubbard model on the 
kagome lattice described above. For any U > 0, the ground states have Stot — Suiax{— M/2) 
and are nondegenerate apart from the trivial spin degeneracy. 

Mielke [35] also extended his results to the situation where the electron density^^ Ai'g/IAI 
is less than 1/3 but close to 1/3. 

In Mielke's work, it was proved for the first time that the Hubbard model with finite 
U can exhibit saturated ferromagnetism. The model is very simple, and the result is very 
important. As far as the author knows, there had been no discussions about the possibility 
of ferromagnetism in the Hubbard model on the kagome lattice. Mielke's work is not only 
mathematically rigorous, but important from physicists' point of view as it opened a new 
way of approaching itinerant electron ferromagnetism. 

Mielke's proof of his main theorem is an elegant induction which makes use of a graph 
theoretic language. The proof is not at all trivial since the problem is intrinsically a many- 
body one. However, there is a very special feature of the model that any ground state of 
the total Hamiltonian H — i^hop + -f^int is at the same time a ground state of each of i^hop 
and i^int- Because of this property, one does not have to face the very difficult problem 
in many-body problems called the "competition between i^hop and i/int-" That one has 
ferromagnetism in this model for any U{> 0) is closely related to this fact. 

Mielke's theorem applies not only to the Hubbard model on the kagome lattice but to 
those on a wide class of lattices called fine graphs. In all of these models, the ground states 
in the corresponding single-electron Schrodinger equation are highly degenerate. There have 
been constructed [36, 37] other examples of Hubbard models in which the corresponding 

^^Thcrc is a minor error in the derivation of critical electron density in Mielke's paper. One should modify 
this part by using the method of [37]. 
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single-electron ground states are highly degenerate, and exhibit saturated ferromagnetism 
for any U > 0. Ferromagnetism in the examples of Mielke and in [36, 37] are now called 
fiat-band ferromagnetism^^. Mielke [38] obtained a necessary and sufficient condition for a 
Hubbard model with highly degenerate single-electron ground states to exhibit saturated 
ferromagnetism. It is interesting that Lieb's ferrimagnetism discussed in Section 5.3 resem- 
bles flat-band ferromagnetism in that the corresponding single-electron spectrum has a bulk 
degeneracy. 

It is needless to say that the models in which single-electron ground states are highly 
degenerate are rather singular. By adding a generic small perturbation to the hopping 
Hamiltonian, the degeneracy is lifted in general, and one gets a nearly flat lowest band 
rather than a completely flat one. A very interesting and important problem is whether 
ferromagnetism remains stable after such a perturbation is added. Of course one does not 
have ferromagnetism for small enough U if the bulk degeneracy in the single-electron ground 
states is lifted. What one expects is that ferromagnetism remains stable when U is sufficiently 
large. (Recall that in the toy model of Section 6.2, we had ferromagnetism for all C/ > only 
for the special choice of parameters t — t'.) There are some indications (from numerical or 
variational calculations) that ferromagnetism is stable under perturbation. As for rigorous 
results, stability of ferromagnetism under single-spin flip is proved in [39, 40] for the model 
obtained by adding arbitrary small perturbation to the Hubbard model of [36, 37]. As for 
a special class of perturbations, the problem of stability of ferromagnetism is completely 
solved as we shall see in the next section. 

6.5 Ferromagnetism in a Non- Singular Hubbard model 

We have seen two theorems which show that certain Hubbard models exhibit saturated fer- 
romagnetism. In Nagaoka's theorem, it is assumed that the system has exactly one hole, 
and has infinitely large Coulomb interaction. In Mielke's theorem and other flat-band ferro- 
magnetism, it is essential that the single-electron ground states have a bulk degeneracy. Is it 
possible to prove the existence of saturated ferromagnetism in a non-singular Hubbard model 
which have finite U and in which the single-electron spectrum is not singular? Recently such 
examples were constructed [41]. 

For simplicity, wc concentrate on the simplest models in one dimension^^. Take the one 
dimensional lattice A = {1,2,..., A^} with A^ sites (where A^ is an even integer), and impose 
a periodic boundary condition by identifying the site A^ + 1 with the site 1. The hopping 

■^°Froin the view point of band structure in the single-electron problem, the bulk degeneracy in the single- 
electron ground states correspond to the lowest band being completely dispersionless (or flat). 

■^"'^ There arc models in higher dimensions [42]. In the original paper [41], the model contains an additional 
parameter A > 0. Here we have set A = V2 to simplify the discussion. The proof of the main theorem in 
[41] is considerably improved in [42]. The condition A > Ac in [41] is replaced by A > 0. 



29 



t t t 




-s -s 



Figure 11: An example of non-singular Hubbard model which exhibits saturated ferromag- 
netism [41]. If we look at three adjacent sites, the lattice structure and the hopping resemble 
that of the toy model of Figure 5. 

matrix is defined by setting tx,x+i = tx+i,x = t' for any x G A, tx^x+2 — ^x+2.a; — t for even x, 
tx,x+2 = tx+2,x = —s for odd X, and tx,y = otherwise. Here t > and s > are independent 
parameters, but the parameter t' is determined as t' = \^{t + s). As can be seen from Figure 
11, the model^^ has two kinds of next nearest neighbor hoppings t and —s, as well as the 
nearest neighbor hopping t'. If we look at an odd site and the two neighboring even sites, 
the model is exactly the same as the toy model we treated in Section 6.2. Roughly speaking, 
this resemblance is the basic origin of ferromagnetism in the present model. We also note 
that because there are next nearest neighbor hoppings, Lieb-Mattis theorem (Theorem 4.1) 
does not apply to the present model. 

The single electron energy eigenvalue in this model can be expressed by using the wave 
number k = 27rn/N {n = 0, ±1, . . . , ± {(A^/4) - 1},N/A) as ei{k) = -2t-2s(l + cos2fc) and 
e2{k) = 2s + 2t{l + cos2fc). There are two bands, and each of them has healthy dispersion. 

As for the Coulomb interaction, we set Ux = U > for any a; G A. We fix the electron 
number as = N/2. In terms of filling factor, this corresponds to the quarter filling. The 
maximum possible value of total spin is S'max = A^/4. 

Unlike in fiat-band ferromagnetism, there is no saturated ferromagnetism when U is suf- 
ficiently small. Theorem 6.1 ensures that the ground state have Stot < •S'max U < 4s. If 
the present system were to show saturated ferromagnetism, it should be in the "nonper- 
turbative" region with sufficiently large U. The following theorem of [41] provides such a 
nonperturbative result. 

Theorem 6.5 (Ferromagnetism in a non-singular Hubbard model) Suppose that the 
two dimensionless parameters t/s and U / s are sujficiently large. Then the ground states have 
Stot = 5'niax(= N/4) and are nondegenerate apart from the trivial spin degeneracy. 

Solvable Hubbard models with U = oo which have similar structure as the present models were found 
by Brandt and Giesekus [43], and were extended in [44, 45, 46, 47]. The conjectured uniqueness of the 
ground state was proved in [48, 47]. The ground state correlation functions in one dimensional models were 
calculated exactly in [48, 49], and insulating behavior was found. ([46] contains an error which is corrected 
in the footnote 6 of [47] . Although I discussed the possibility of superconductivity in these model in [46] , I 
am not very optimistic about this conjecture at present.) 
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The theorem is vahd, for example, when t/ s > A.h \l U / s — 50, and t/s > 2.6 if 
U/s — 100. The ferromagnetic ground state can be constructed in exactly the same manner 
as $1 in the previous section. 

Although the model is rather artificial, this is the first rigorous example of saturated 
ferromagnetism in a non-singular Hubbard model on which we have to overcome the com- 
petition between ifint and i^hop- In ^ class of similar models, it is also proved that low-lying 
excitation above the ground state has a normal dispersion relation of a spin-wave excitation 
[41, 40]. Starting from a Hubbard model model of itinerant electrons, the existence of a 
"healthy" ferromagnetism is established rigorously. 

If we set s = in the present model, the ground states of the single-electron Schrodinger 
equation become A^/2-fold degenerate. In this case, the model exhibits saturated ferromag- 
netism (fiat-band ferromagnetism) for any U > 0. Theorem 6.5 for s ^ can be regarded as 
a solution to the problem about stability of flat-band ferromagnetism against perturbation 
to the hopping Hamiltonian. 

The basic strategy in the proof of Theorem 6.5 is first to estabhsh the existence of 
saturated ferromagnetism in a Hubbard model on a chain with five sites, and then "connect" 
together these local ferromagnetism to get ferromagnetism in the whole system. Generally 
speaking, this is a crazy idea. In a quantum mechanical system, especially in a system with 
healthy dispersion (like the present one), electrons have strong tendency to extend in a large 
region and reduce kinetic energy. To confine electrons in a finite region usually costs extra 
energy. To obtain exact information about large system from smaller system seems to be 
impossible. The reasons that such a strategy works in the present model arc twofold. One is 
the special construction of the model. The other is that wc described electron states using a 
language which takes into account both the particle-like character of electrons and the band 
structure of the model. The latter is a natural strategy to deal with the Hubbard model in 
which wave-particle dualism generate interesting physics. 

It is believed that the ground states of the present Hubbard model describes an insulator. 
When the electron number is less than N/2, we expect that the present model exhibits 
metallic ferromagnetism in which the same set of electrons participate in conduction as well 
as magnetism^^. For the moment, we still do not know of any useful ideas in proving this 
fascinating conjecture. 

I wish to thank Kenn Kubo and Balint Toth for useful comments on the early version of 
the present article. 



By using a heuristic perturbation theory based on the Wannier functions, the low energy effective theory 
of the present Hubbard model is shown to be the ferromagnetic t-J model. Moreover, by considering models 
close to the flat-band model, one can make \t/J\ arbitrarily small. This observation gives a strong support 
to the above conjecture. 
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